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Abstract. We analyze optical and near infrared spectra of intermediate and high resolution of the eclipsing 
symbiotic system AR Pavonis. We have obtained the radial velocity curves for the red and the hot component 
from the M-giant absorption lines and from the wings of Ha, H/3 and He II A4686 emission profiles, respectively. 
From the orbital elements we have derived the masses, M g = 2.5M© and Mh = I.OMq, for the red giant and 
the hot component, respectively. We also present and discuss radial velocity patterns in the blue cF absorption 
spectrum as well as various emission lines. In particular, we confirm that the blue absorption lines are associated 
with the hot component. The radial velocity curve of the blue absorption system, however, does not track the 
hot companion's orbital motion in a straightforward way, and its departures from an expected circular orbit are 
particularly strong when the hot component is active. We suggest that the cF-type absorption system is formed 
in material streaming from the giant presumably in a region where the stream encounters an accretion disk or an 
extended envelope around the hot component. The broad emission wings originate from the inner accretion disk or 
the envelope around the hot star. We also suggest that the central absorption in H profiles is formed in a neutral 
portion of the cool giant's wind which is strongly concentrated towards the orbital plane. The nebula in AR Pav 
seems to be bounded by significant amount of neutral material in the orbital plane. The forbidden emission lines 
are probably formed in low density ionized regions extended in polar directions and/or the wind- wind interaction 
zone. 
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(1937) as an eclipsing binary with a period of 605 days. 
The eclipsed object is highly variable in both brightness 



1. Introduction citation in the eclipsed region increasing towards the cen- 

tral primary light source. TH74 derived masses of order 
The s ymbiotic star AR Pavonis was discovered by Mayall 2 .5M Q and 1.2 M , for the primary and secondary, re- 
spectively. Their estimate was based on analysis of the 
eclipse shape, and the primary mass function derived from 
radial velocities of almost exclusively Hei emission lines. 
In particular, they did not measure radial velocities for 
the cool giant. TH74 also found periodic changes in radial 
velocities of forbidden lines and cF-type absorption lines, 
which they attributed to a stream and a region where 
the stream falls on a ring /disk surrounding the primary. 
Kenyon & Webbink ( 1984 ) revised this model involving an 
accretion disk around a main sequence star as the primary 
light source. 



and size (Andrews 1974). Thackeray & Hutchings (1974, 
hereafter TH74) based on extensive analysis of spectro- 
scopic data spanning several orbital cycles proposed a bi- 
nary model for AR Pav with the M3 III secondary fill- 
ing its Roche lobe. They found that the eclipsed primary 
shows an emission O-type spectrum and density and ex- 
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Bruch et al. (1994) observed dramatic variations in the 
visual light curve between eclipses which are not easily ex- 
plained by the TH74 model. They suggested variations in 
the mass transfer rate from the red giant to the hot com- 
ponent as the most satisfying explanation. Skopal et al. 
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Table 1. Log of the spectroscopic observations of AR Pav 
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( 2000a ) noted that the AR Pav light curve has a charac- 
teristic shape similar to those of dwarf novae in quiescent 
phase; whereas during the active phases the profile of the 
light curve changes considerably. 

The first orbital solution for the cool giant has been 

hereafter S01). 



(2001 



recently obtained by Schild et al. 
They also measured the rotation velocity of the giant, and, 
assuming co-rotation, its radius. Based on this radius com- 
bined with the known spectral type of the giant and the 
corresponding effective temperature S01 estimated the gi- 
ant's luminosity and distance to AR Pav, as well as from 
the position in the HR diagram the giant's mass. Their 
binary mass function then resulted in masses of 2 Mq and 
0.75 Mq for the cool giant and its hot companion, respec- 
tively. 

In this paper we present new radial velocity data for 
the cool giant and the emission lines collected between 
1990 and 2000. From this we derive new spectroscopic 
orbits for both components, and discuss the velocity pat- 
terns in various emission lines. 

2. Observations 

Spectroscopic observations were performed with the 
2.15 m telescope of Complejo Astronomico el Leoncito, 



CASLEO (San Juan, Argentina). During 1990 and 1991 
medium-low resolution spectra were taken with a Boiler & 
Chivens Cassegrain spectrograph using a photon-counting 
Reticon, called Z-Machine (A/AA=4100 and 2700 in the 
filters blue and red respectively), and a Thompson CCD 
of 384x576 pixels (A/AA=2600) in 1992. Since 1995, 
high resolution spectra were obtained with a REOSC 
echelle spectrograph using a Tek CCD 1024x1024 pixels 
(A/AA=15000), except in May 1998 when a spectrograph 
Boiler & Chivens was used. Table |l| shows a log of the 
observations. 



Z M: Reticon Z-Machine 

T CCD: Thompson CCD (384x576 pixels) 

R CD: REOSC echelle spectrograph in cross-dispersion 

BC: Boiler & Chivens spectrograph 

1 images taken with bin factor=2 

2 image taken with bin factor=4 



Barba et al. (1992) describe details of the Z-Machine, 
acquiring and reducing data procedures. The CCD data 
were reduced with IRAFQ packages, CCDRED and 
ECHELLE and all the spectra were measured using the 
SPLOT task within IRAF. 

To obtain the flux calibration, standard stars from 
Stone & Baldwin (p83|) and Baldwin & Stone rtl98l 



were observed each night. A comparison of the spectra 
of the standards suggests that the flux calibration errors 
are about 15 per cent for the Z-Machine and 20 per cent 
in the central part of each order for the REOSC echelle 
images. 

In the determination of the hot component orbit, we 
have included the Ha and He IIA4686 profiles measure- 
ments from a spectrum taken by Van Winckel ( 1993 ), in 
July 1988 and corresponding to phase <fi — 0.72. 

3. Analysis and discussion 

3.1. Emission line profiles 

The left panel of Fig. |] shows the Ha emission line profiles 
of AR Pav ordered with their orbital phase, although they 
are collected from more than one cycle. Ha is strongest in 
phase 0.54, when the hot component is in front. The pro- 
files are double-peaked with broad, FW <; 1600 km s , 
wings. The red peak is stronger than the blue one in all 
phases. This can not be explained by motions in an op- 
tically thin medium. In such situation the profile of the 
second quadrature would be reversed relative to the first 
quadrature. The profiles are apparently affected by self- 
absorption. The largest blue to red peak intensity ra- 
tio is observed around phases 0.13 and 0.27, respectively, 
whereas the smallest value of this ratio corresponds to 
phase 0.77. The double-peaked profiles have an average 
separation of 148 ± 3kms _1 . The average radial velocity 
of the central absorption is —96 ± lkms -1 . In eclipse, 
Ha profiles are considerably narrower than out of eclipse 
(Fig. H|) . This suggest that the broad Ha wings are formed 
near to the hot component. 

The shape of H/3 profile is similar to that of Ha, but 
the blue to red peak intensity ratio is smaller than the 
same ratio for Ha, except for phases 0.98 to 0.13, when 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which is operated by the association of 
Universities for Research in Astronomy, INC., under contract 
to the National Science Foundation 
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Fig. 1. Emission line profiles in AR Pav. The continuum was subtracted and each profile was normalized to maximum 
intensity. The profiles are shifted vertically by 0.5 (Fell) and 0.4 (other lines), respectively, for better display. Ha 
profile corresponding to August 95 is saturated. In the right panel, solid line shows Fe II A5316 and dashed line shows 
Fe II A5363. Only high resolution spectra are shown. 



the relative intensities are similar. The wings are also nar- 
rower during eclipse. The double-peaked profiles have a 
average separation of 117 ± 3kms~ 1 , and the average ra- 
dial velocity of the central absorption is —90 ± 2 kins" 1 . 

We have H7 profiles only in phases near to both con- 
junctions. The intensity of the blue peak relative to the red 
one is smaller than for H/3, except during eclipse. It seems 
to behave in the same way as Ha and H/3. The average 
radial velocity of the central absorption is — 87±3kms™ 1 . 
The absorption components of the Balmer lines seem to 
show a positive progression from Ha to H7. Similar posi- 
tive progression in radial velocities of the central absorp- 



tion component of Balmer emission lines is often observed 
in symbiotic stars. Since the absorption lines have slower 
decrement than the emission lines the positive progression 
may be due to the decreasing effect of the emission border. 

The second panel of Fig. |l| shows the He 11 A 4686 emis- 
sion lines profiles. The line profile is double-peaked during 
eclipse and in phase 0.85. During eclipse the blue peak is 
slightly stronger than the red one. In phase 0.85 the red 
peak is stronger than the blue one. In the rest of the or- 
bital phases they present asymmetric profiles. The asym- 
metry is reversed in the first quadrature relative to second 
one. However, in phase 0.13 the shape of the asymmetry 
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Fig. 2. The wings of Ha at selected orbital phases. 



is almost identical to phase 0.77. The largest asymmetry 
occurs in phase 0.54. 

Fig. |] also shows the [O in] A 5007 emission line re- 
gion. The permitted lines of Hei A 5016 and Fen A 5018 
are narrow and more or less symmetric, with FW ~ 
100 - 150 km s" 1 . The forbidden [Om] lines display FW 
~ 200 km s" 1 and may have double peak structure with 
peak-to-peak separation of Av ~ 40kms _1 , and the blue 
to red peak intensity ratio of ~ 1. This structure is visible 
in our profiles at phases 0.54 and 0.89 for which we esti- 
mate Aw ~ 40 and 50kms _1 , respectively (Fig. [j]), as well 
as in the profiles presented by van Winckel et al. ( ^993| ) 
and S01. 



In the high resolution spectra, also some Fell emis- 
sion lines appear as double-peaked and/or asymmetric at 
most orbital phases (the right panel of Fig. |l|). In partic- 
ular, the profiles are clearly double-peaked near the first 
quadrature, and asymmetric (with a very weak central ab- 
sorption occasionally present) at the second quadrature. 
The profiles seem symmetric and the central reversal is 
missing during the eclipse whereas at <fi ~ 0.54 a very 
faint blue peak (hardly discernible from the noise) can 
be present. The double-peaked Fell profiles have separa- 
tions between 47 to 65 kms -1 . In the spectra of June 
1995, the He I emission lines AA 4922, 5876 and 6678 show 
also a double-peaked profile with an mean separation of 
34 kms^ 1 , whereas the Hei A 5016 an d the Fen A 5018 pro- 
files published by van Winckel et al. ( 1993 ) show a reversal 
at ~ —93 and ~ — 86kms~ 1 , respectively. It is remarkable 
that in both conjunctions the intensity of Fen A 4923 and 
A 5018 emission lines are comparable with the intensities 
of He I A 4921 and A 5015, respectively. In other phases the 
ratio Fe/He is very small. 



3.2. Radial velocities 

To obtain the radial velocities of the red giant and the hot 
component, we have measured spectral features associated 
with each one. 

For the cool component, we have used M-type absorp- 
tion lines corresponding to Fe I, Ti I, Ni I, Si I, O I, Zr I, Co I, 
Vi, Mgi and Gdli. At some phases we have also identi- 
fied and measured the blue cF-type absorption lines cor- 
responding to Cm, Fen, Tin and Yn which are believed 
to be linked to the hot companion (TH74; S01). The indi- 
vidual radial velocities were obtained by a Gaussian fit of 
the line profile, and a mean value was calculated for each 
spectrum. The resulting mean heliocentric velocities to- 
gether with their standard errors and the number of lines 
used are given in Table ^. 

In addition, for the study of the hot component we 
have determined the radial velocities from the wings of 
Ha, H/3 and He II A 4686 profiles. They would reflect the 
motion of the hot component if they are formed in the 
inner region of the accretion disk or near to the hot com- 
ponent. In this way we have used a method outlined by 
Schneider & Young (1980) and Shatter ( 1983| ), which was 
very successful in studies of cataclysmic binaries. This 
method consists of convolving the data (S) with two iden- 
tical Gaussian bandpasses (G) whose centers have a sep- 
aration of 2b. The wavelength, A, of an emission line in a 
spectrum S(A) is given by solving the equation 



5(A)G(A - A)dA = 



where 

G{x) = exp[-(x - b) 2 /2a 2 } - exp[-(a; + b) 2 /2a 2 ]. 

The choice of the parameters b and a depend on the char- 
acteristic of the spectra being analyzed, namely the emis- 
sion line width and the si gnal-to-noise ratio of the data. 
We refer to Shatter ( 1983 ) for a more detailed discussion 
of the influence of the choice of b and a on the obtained 
orbital solution. In any case they found that varying these 
parameters has no significant effect on the resulting value 
of the system velocity, 7, whereas the semi-amplitude of 
the radial velocity, K, significantly depends on the choice 
of b. The resulting velocities are given in Table |[ 

In order to reduce the noise in the extreme line wings, 
we have added all spectra taken in the same observing 
run. Fig. [l] and Fig. || show the Ha emission line profiles 
and wings, respectively, of AR Pav ordered with their or- 
bital phase. We have experimented with a range of val- 
ues for &, and a fixed at 7 A, for which we have obtained 
the circular orbit solutions (see below). Fig. || illustrates 
the dependence of the semi-amplitude, K, and its stan- 
dard error, ctk on b in the case of Ha. It is noticeable 
that ck increases sharply for b larger than ~ 10 A. We 
can attribute these large standard errors to the velocity 
measurements dominated by the noise in the continuum 
rather than by the extreme high velocity wings of the line 
profile. We have therefore adopted b = 10 A(460 kms -1 ) 
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Fig. 3. The orbital semi-amplitude K and its standard 
error ok for the wings of Ha (dots) and H/3 (open circles) , 
respectively, as a function of the parameter b. The best 
estimation of K is determined by the value of b where ok 
begins to sharply increase 



as the best value in the case of Ha. A similar analysis led 
us to choose b — 7 A (430 km s" 1 ) and a = 5 A for H/3, 
and b = 5 A (320 km s" 1 ) and a = 3 A for Hen A 4686, 
respectively. 

We have also obtained the radial velocities of He I AA 
4009, 4026, 4121, 4144, 4388, 4471, 4713, 4922, 5016, 5048, 
5876, 6678, 7065, 7281; [Om] AA4363, 4959, 5007; [Nil] 
A 6584 and occasionally also A 5755 emission lines as well 
as numerous Fell emission lines by a Gaussian fit. Faint 
[Fen] though present at some phases have not been mea- 
sured. Table |2| shows the average heliocentric radial ve- 
locities, their errors and number of lines used for each 
spectrum. The singlet and triplet series of Hei lines are 
shown separately. In the case of double-peaked lines (see 
above) of Fe II the center velocity were used to derive the 
values in Table ||. 

3.3. Spectroscopic orbits 

All radial velocities in Table |^ except the central absorp- 
tion of H I vary with the 605-day period. Fig. ^ shows our 
radial velocity data corresponding to the M giant absorp- 
tion lines, the H I and He II emission wings, and He I, Fe II, 
[Om] and [Nil] emission lines, respectively, phas ed wi th 
the photometric ephemeris given by Bruch et al. ( 1994 ) 
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Fig. 4. Radial velocity data and circular orbital solution 
for AR Pav. Closed symbols represent our data, and tri- 
angles represent the data of S01 (M giant) and TH74 (cF 
absorption and emission lines). The solid line repeats the 
orbit of the M giant and the dotted line - the hot com- 
ponent solution, respectively (a) M giant absorption, (b) 
Hot component (HC): filled circles correspond to the Hi 
wings, filled squares correspond to Hen A4686. (c) cF ab- 
sorption lines, (d) Hei emission lines, (e) Fell emission 
lines, (f ) [O in] emission lines, (g) [N n] emission lines. 



Our radial velocity curve for the M giant is practically 
indistinguishable from that given by the data from S01 
(and also plotted in Fig. ^), and it strongly suggests a 
circular orbit. It is also obvious from Fig. Q that the cF 
absorption system, H I and He n emission wings and the 
permitted (He I and Fe n) emission lines are in antiphase 
with the M giant absorptions. 

In order to obtain the orbital parameters of AR Pav 
we have assumed a circular orbit as did S01, and adopted 
the very precise photometric ephemeris of Bruch et al. 



JDMin = 2 420 331.3 + 604.5 E. 



(1994). The semi-amplitude of the respective radial veloc- 
ity curve has been then fitted by a least squares method. 
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Table 2. Radial velocities of the giant absorption lines, the blue cF-type absorption lines, permitted and forbidden 
emission lines, emission wings of Hi and Hen lines*, and the central absorption component of Ha (HA) in AR Pav 
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-71+1 (3) 


-89(1) 








-100 


51629/0.774 


-80+1 (26) 


-43+1 (17) 


-44+1 (29) 


-47+1 (5) 


-44+2 (4) 


-68+0 (2) 


-79+3 (2) 


-34 


-41 


-33 


-93 


51631/0.778 


-79+3 (21) 


-42+1 (15) 


-47+2 (27) 


-48+1 (5) 


-44+3 (4) 


-69+1 (2) 


-85(1) 








-94 


51755/0.982 


-70+3 (14) 




-60+1 (33) 


-63+2 (6) 


-60+2 (5) 


-63+1 (3) 


-82+4 (2) 


-56 


-61 


-59 


-93 


51757/0.986 


-68+1 (66) 




-62+2 (13) 


-59+1 (2) 


-61+2 (2) 




-82+3 (2) 
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★ Sum of the s pectra taken in the same observing run. 
f Van Winckel( |l993| ) 
s: saturated emission 



In some cases we have also calculated elliptical orbital so- 
lutions with the orbital period forced to the photometric 
period of 604.5 d. Table || lists the resulting orbital so- 
lutions. The symbols have their usual meaning: 7 is the 
system velocity, K is the orbital semi-amplitude, e is the 
eccentricity, uj is the longitude of periastron, To is the 
time of periastron passage, f(M) is the mass function, 
Asini is the fractional semi-major axis in astronomical 
units (AU), and AT is the time difference between orbital 
conjunction and photometric eclipse. For the M giant, we 
have obtained both solutions using only our radial veloci- 
ties as well as for our radial velocities combined with the 
S01 data. Although an elliptical orbit fits the measured 
velocities slightly better than a circular one, the obtained 
eccentricity e = 0.12 + 0.04 is only a 3a result, and ob- 
servations covering several more orbital periods would be 
necessary to confirm the reality of such an orbit. Following 
the arguments given by S01 we believe that a circular or- 
bit provides a better description of the data, and such an 
orbit is also in agreement with the tidal theory. 

Both the broad emission line wings, and the cF-type 
absorption lines show the highest (and practically the 
same) amplitude and the mean velocity almost identical 
with the red giant systemic velocity. The solutions for 
the broad emission wings of Hi and He 11, respectively, 



agree within their respective errors which suggests that 
both they are formed in the same region near the hot 
component. An A- or F-type absorption system with ra- 
dial velocity varying in antiphase with the M giant ab- 
sorption lines, similar to that identified in AR Pav, was 
seen during outbursts of two other symbiotic systems AX 
Per and BX Mon, and associated with the hot compo- 
nent (Mikolajewska & Kenyon 1992a, Dumm et al. 1998). 



This combined with the fact that the wings narrow dur- 
ing the hot component eclipse strongly suggests that both 
the broad emission wings and the blue (cF) absorption 
system trace the orbit of the hot component. The only 
problem is that any orbital solution for the cF absorption 
system leads to significant eccentricity (Table [| TH74; 
S01) whereas the red giant orbit requires a circular (or 
nearly circular) orbit for the hot companion. Our orbital 
solution for the broad emission line wings, on the other 
hand, is consistent with a circular orbit. 

To analyze the cF absorption system we have com- 
bined our measurements with the TH74 data. Since our 
measurements are a factor of 2-4 more accurate then those 
of TH74, we have applied a weight 4 to our data, and 
weights 1-3 to the TH74 data (same as in their Table 1), 
respectively. We have re-calculated the orbital solution for 
the combined data and found practically the same eccen- 
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Component 


7 [km s J 


i\ [Km s 


e 






/ \ M ) I m q\ 


J\ sin I [A. U J 


AT( 2 ) 


M abs (our) 


-68.3 ±0.8 


11.4 ± 1.2 


775T 

qW 












M (our+SOl) 


-68.4 ±0.3 


10.8 ±0.5 


q(3) 






0.079 ±0.011 


0.60 ±0.03 


(3) 


M (our±S01) (4) 


-68.2 ±0.2 


10.9 ±0.4 


0.12 ±0.04 


351 ± 12 


50693 ± 20 


0.079 ± 0.009 


0.60 ± 0.02 


12 


cF (our+TH74) (5) 


-66.5 ± 1.3 


27.6 ±2.9 


0.45 ± 0.07 


273 ± 11 


50899 ± 21 


0.94 ± 0.40 


1.37 ±0.19 


34 


cF (our±TH74) (6) 


-67.4 ± 1.1 


22.8 ± 1.7 


0.21 ±0.07 


307 ± 20 


50936 ± 31 


0.69 ±0.18 


1.23 ±0.11 


4 


Wings (Hi) 


-65.2 ± 1.9 


26.3±2.7 


(3) 






1.14 ±0.32 


1.48 ±0.15 


(3) 


Wings (Hi) (4) 


-64.8 ± 1.4 


25.8 ± 1.8 


0.13 ±0.07 


68 ±29 


51726 ± 45 


1.05 ± 0.26 


1.42 ±0.12 


-7 


^A/incc TT^ 
VV lllgjO ^llc 11 j 


— 70 A + 9 7 




(3) 






1 ^ 1 4- fl 70 

l.Jl ZC U. 1 u 


i fin 4- n 99 


(3) 


Wings (Hi+Heii) 


-67.1 ± 1.6 


27 .4 ±2.3 


(3) 






1.29 ±0.35 


1.52 ±0.13 


Q(3) 


He i 


-63.3 ±1.8 


15.0 ±2.8 


Q(3) 










Q(3) 


He i (s) 


-64.2 ± 1.8 


15.6 ±2.8 


Q(8) 










Q(3) 


He I (t) 


-62.2 ± 2.0 


14.1 ±3.1 


Q(3) 










Q(3) 


Fen 


-63.6 ±0.9 


17.8 ± 1.4 


Q(8) 










Q(3) 



(1) the time of periastron passage; 

(2) AT = Tconj — Tecl ipse j 

(3) assumed; 

(4) unconstrained solution; 

(5) unconstrained, weighted solution; 

(6) same as (5) for quiescent data only (see text). 



tricity, e = 0.45 ± 0.07, as e ~ 0.45 and 0.42 given by 
TH74 and S01, respectively, and 7 = — 66.5kms _1 , closer 
to the red giant systemic velocity than earlier solutions. 
The spectroscopic conjunction occurs 34 days after the 
photometric eclipse. S01 suggested that the cF absorption 
system only traces the hot component's motion whilst it 
is in front of the giant. At other orbital phases, the cF 
absorption lines are affected by additional absorption due 
to this wind and/or a flow of material towards the hot 
companion, and the resulting radial velocities mimics an 
elliptical orbit. S01 have also demonstrated that a circu- 
lar orbit that is almost exactly in antiphase with the red 
giant results if only the radial velocity measurements be- 
tween phase 0.25 and 0.75 ( the hot component in front) 
are employed. Unfortunately, our new measurements, in 
particular 2 points almost coincident at phase ~ 0.5 do 
not fit this restricted orbital solution of SOI. Instead we 
have found that the cF absorptions may be affected by 
the hot component activity. This is illustrated in Fig. 0. 
In particular, all data points that force an elliptical orbit - 
the TH74 points between phase ~ and 0.2 as well as our 
two points near phase ~ 0.5 - were obtained when AR Pav 
was very active (TH74; Skopal et al. pOOObl) . It is interest- 
ing that an orbital solution omitting the TH74 data from 
the 1954 outburst as well as our measurements obtained 
prior to JD ~ 2 451000 (when according to Skopal et al. 



-30 



(2000b) AR Pav was very active) results in a much lower 
eccentricity, e = 0.21 ±0.07 and the spectroscopic conjunc- 
tion at the time of the photometric eclipse. Both these fits, 
however, match the data rather poorly, and the residuals 
for some points are larger then their measurement accu- 
racy. This illustrates difficulties with interpretation of the 
radial velocity curve of the blue absorption system. 



cF-abs 




9 our - active 
■ our - quiet 
n TH74 - active 
A TH74 - quiet 



0.5 1 

Phase 

Fig. 5. The radial velocity curve for the cF shell com- 
ponent with our circular orbits of the red giant (solid 
curve) and the broad emission wings (dashed curve), re- 
spectively. The horizontal line represents the system ve- 
locity, 7 = — 68.4km.s~ 1 . Open symbols (triangles and 
stars) represent the TH74 data, with sizes reflecting their 
weights. Filled symbols (squares and dots) correspond to 
this work. We use different symbols for data obtained dur- 
ing activity and quiescence, respectively. The dotted curve 
gives an elliptical fit to all data points, the dot-dashed 
curve - to quiescent data only. 



We note here that similar complications with the blue 
absorption system occurred for a few other active symbi- 
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otic systems. In the case of symbiotic recurrent nova RS 
Oph, Dobrzycka & Kenyon ( |1994[ ) found the A-type ab- 
sorption features ~ 0.4 P out of phase with the orbit of 
the giant. Mikolajewska & Kenyon ( |1996 ) failed to derive 
any radial velocity curve and orbital solution for the blue 
absorption system in Z And, whereas in CI Cyg the radial 
velocities of the F-type absorption system suggest forma- 
tion region in material streaming from the giant near the 
hot component (Mikolajewska & Mikolajewski 1988). 

Such problems are also common in the case of W Ser 
stars (or hyperactive Algols), in which the accreting pri- 
mary is embedded in an optically thick shell with an A/F- 
type spectrum, and strong emission from circumstcllar 
plasma is also present. We recall here that the radial ve- 
locity curves for the shell spectra in W Ser stars show of- 
ten significant departure from that expected for a circular 
orbit required by the orbital solut ions for their secondary 
components (Andersen et al. 1988 ) . It has been also found 
that such significant departures can be expected if the 
shell lines used to measure the radial velocity curve are 
forme d in t he outer 10 — 20% of the disk (e.g. Andersen 
et al. 1988). In particular, the distorted radial velocity 
curve of the shell in SX Cas (Fig. 1 of Andersen et al. 



1988f ) is very similar to the radial velocity curve of the 



blue absorption system in AR Pav. 

The possible correlation between AR Pav's activity 
and the departures of the cF absorption velocities from the 
circular orbit suggests that this absorption system may be 
also affected by material streaming towards the hot com- 
ponent presumably in a region where the stream encoun- 
ters an accretion disk or an extended envelope around the 
hot component. The outer disk and the stream impact re- 
gion can be highly unstable and asymmetric. For example, 
the trailing side of the disk, where the gas stream adds to 
the disk, can be brighter whereas the leading side is more 
extended. Such asymmetry may, at least qualitatively, ac- 
count for the deviations of the cF absorption velocity curve 
from the sinusoid expected for a circular orbit as well as 
for the behaviour of the optical light curves of AR Pav, 
which show relatively stable profile of eclipses during qui- 
escence in contrast to very unstable profiles and changes 
in th e position of the m inima durin g activ ity (Andrews 
1974 Bruch et al. |1994| ; Skopal et al. ^000a| ). 



Unfortunately, a detailed and proper treatment of the 
cF shell radial velocity curve is beyond the scope of our 
paper. This would require knowledge of the temperature 
and density profile in the shell/disk which is not avail- 
able. The situation is further complicated by significant 
cycle-to-cycle variations of the cF spectrum. In particu- 
lar, spectroscopic observations a few orbital cycles and of 
good phase coverage in each cycle are necessary to dis- 
tinguish between the cycle-to-cycle changes of the radial 
velocity curve and distortions due to other effects (e.g. 
contribution from stream, formation in outer disk, etc.). 

The solutions for the broad emission line wings are free 
from such complications and they are not affected by the 
system activity. We thus believe that the broad wings are 
formed very close to the hot star (in the inner disk or in 



a hot wind) and they directly trace the hot component's 
orbit. 

Combing the semi-amplitudes of the M giant and 
the Hi+Hcii emission line wings for the circular orbit 
(Table ||) gives a mass ratio q — 2.5 ± 0.3, and compo- 
nent masses of M g sin 3 i — 2.5 ± 0.6 Mq and Mh sin 3 i = 
1.0±0.2 Mq, with the errors given by the errors of the K^ 
and K g values, respectively. These values are somewhat 
larger than M g sin 3 i = 1.9 Mq and M h sin 3 i = 0.75 Mq 
derived by S01 mostly because our semi-amplitude for 
the M giant, K g = 10.8 km s -1 , is larger than the value 
K g = 9.6 km s -1 calculated by S01. We, however, note 
that our radial velocity curve (Figure ^) has better cover- 
age at both quadratures. 

The radial velocities of the central absorption compo- 
nent of Hi Balmer lines do not follow the orbital phase, 
and are blueshifted by ~ 30kms _1 relative to the sys- 
tem velocity. The absorption is presumably formed in the 
neutral portion of the cool component's wind. Moreover, 
it remains strong and practically constant over whole or- 
bital cycle, which indicates that the nebula in AR Pav is 
probably bounded on all sides by significant amount of 
neutral material at least near the orbital plane. It is in- 
teresting that most symbiotic systems with double-peaked 
Ha profiles have the main reversal on the blue side of the 
line center, strongly supporting the idea that this struc- 
ture arises from self-absorption in the cool giant's wind 



(Ivison et al. |1994J). 

For a weak, unsaturated absorption the relationship 
between column density, N[, and the measured equivalent 
width, EW, is 



Ar 1.13 x 10 20 EW 2 

A /ij 

where /y is the oscillator strength of the transition i —> j, 
and wavelength A and EW are given in A. From the mea- 
sured EW(Ba) = 1.17±0.03A, we estimate the lower 
limit for N(Rl,n= 2) £ 4.4 x 10 12 , and from EW(H.(3) = 
0.66 ±0.06 A we obtain N(R I , n = 2) > 3 x 10 13 cm" 2 , re- 
spectively. To estimate the total hydrogen column density, 
AT(H), the concentration of hydrogen atoms at the second 
excitation level has to be adopted which requires detailed 
knowledge on physical conditions in the gas. Assuming 
local thermodynamic equilibrium (LTE) and a tempera- 
ture, 6000 K, we obtain (from the Saha-Boltzman rela- 
tion) AT (H i,n = 2)/AT(H) k 10~ 8 , and the total hydrogen 
column density AT(H) <; a few x 10 21 cm -2 . Despite of all 
uncertainties involved in this estimate, the resulting hy- 
drogen column density is consistent with the behavior of 
the shortwavelength IUE spectra of AR Pav. These spec- 
tra show characteristic attenuation of the continuum flux 
for A 1500 A at any phase - including SWP 5828 and 
SWP 13956 taken at orbital phases 0.27 and 0.37, respec- 
tively - which can be accounted for by Rayleigh scattering 
in the neutral material surrounding the binary system. It 
has been shown that to produce a noticeable effect col- 
umn densities in the range AT (Hi) ~ 10 20 



10 24 cm 2 are 



required (e.g. Schmid 1997). We thus conclude that both 
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the strength of the central absorption in H I Balmer lines 
and the behavior of the IUE SWP spectra are consistent 
with significant amounts of neutral material bounding the 
nebula in the orbital plane. 

The radial velocity curves for Hei and Fen emission 
lines are in antiphase with the cool giant, which suggests a 
possible association of these lines with the hot component. 
However, their semi-amplitude, ~ 15 and ~ 18kms _1 for 
He I and Fen, respectively, are a factor of ~ 2 lower than 
that of the hot component, and their mean velocities are 
redshifted by ~ 5kms _1 with respect to the system ve- 
locity. Such discrepant amplitudes can be explained by a 
contribution from material between the two stars and in 
the neighborhood of the red giant which can both reduce 
the amplitude and shift the mean velocity (sec also discus- 
sion in SOI). We have not found any significant difference 
between the K values as well as 7 velocity derived for the 
singlet and triplet Hei emission series, respectively. There 
is also no difference between our radial velocity curve and 
the data from TH74, which means that the Hei formation 
region was the same at the two epoches. 

The radial velocities of the forbidden lines show a more 
complicated pattern (Fig. [|, bottom). The [Oni] emission 
lines do not trace any of the two stellar components. Our 
radial velocity data show a deep minimum around the 
phase ~ 0.7 and another less marked minimum around 
the phase ~ 0.2, and the mean velocity, — 68kms , very 
close to our red giant 7 velocity. A small amplitude peri- 
odic RV changes are also present in TH74 data. SOI sug- 
gested that the variability may be due to variable relative 
size of the blue and red peaks causing the line center of 
unresolved [O ill] lines (such as [O in] A 4363 in TH74) to 
move bluewards and redwards. However, we have found 
the same RV values and changes by measuring the line 
centers and the wings, respectively, which in our opinion 
suggests that the changes are due to true motions. 

Finally, the [N 11] emission lines show a small periodic 
changes roughly in phase with the hot component and the 
mean velocity blueshifted by ~ 16kms . 

3.4. The M giant, the hot companion and the 
geometry of AR Pav system 

AR Pav is well-established eclipsing system. The eclipses 
are however not total as indicated by relatively blue 
colours during mid-eclipse, B — V ~ 0.6, much different 
from B — V ~ 1.5 expected for an M5 giant, and variable 
U BV magnitudes during minimum which increase with 
the magnitudes outside eclipse (Andrews 1974; Hutchings 
et al. 1983). Similarly, the UV continuum flux measured 
from IUE spectra (cf. Fig. 5 of SOI) shows a significant re- 
duction near = 0, the continuum is, however, not fully 
eclipsed, and the shape of the light curve is asymmetric 
and complex. SOI estimated a limit i il 79° from the ob- 
served exit from eclipse for a point-like source and using 
the red giant radius derived from its rotational velocity. 



An independent estimate of the orbital inclination can 
be made by comparing the apparent red giant radius in- 
ferred from eclipses, i? g , e ci, with its mean tidal radius, Rt, 
defined as the radius of a sphere with a volume equal to 
the Roche lobe volume. The inclination is 



cos? < 



\l {Rt/Af - (i? g , ecl M)2. 



Skopal et al. fl2000a| ) obtained R^ cc \/A = 0.30 ± 0.020, 
while R t /A = 0.46 ± 0.01 for q = 2.5 ± 0.3 (Paczynski 



1971 ). If the giant actually fills its tidal lobe, the orbital in- 



clination is i = 70 ±2 , so M g ~ 3 M and M h ~ 1.2 Mq. 
Otherwise these estimates give a strong lower limit for i 
and upper limits for the binary component masses. 

Our orbital solution yields a binary separation of 
Asini = 457±35 Rq. For i = 90°, the giant radius derived 
from eclipses is then 0.3 A — 137±20Rq, very close to the 
value 130 ± 25 Rq derived by S01 from rotational veloc- 
ity. The red giant radius extends to ^ 65 % of the mean 
Roche lobe radius^] with the lower limit set by i = 90°. 
Based on similar estimate for the Roche lobe filling factor 
S01 concluded that the giant in AR Pav does not fill its 
tidal lobe. In our opinion, however, this result does not 
settle the problem. 

Recently, Orosz & Hauschildt ( |2000 ) have shown that 
rotational broadening kernels for Roche lobe filling (or 
nearly filling) giants can be significantly different from 
analytic kernels due to a combination of the nonspheri- 
cal shape of the giant and the radical departure from a 
simple limb darkening law. As a result, geometrical in- 
formation inferred from v sin i measurements of cool gi- 
ants in binary systems, and in particular in symbiotic 
stars, are likely biased and must be treated with cau- 
tion. In particular, among symbiotic systems with mea- 
sured vsm i and accurate orbits, three of them, T CrB, 
CI Cyg and BF Cyg, seem to contain tidally distorted gi- 
ants as indicated by their red/near- IR light curves with 
evident ellipsoidal changes (Belczynski & Mikolajewska 
19981 Mikolajewska |200l[ Mikolajewska et al. |2002|). In 



CI Cyg, the rotational velocity of the giant is consistent 
with a synchronously rotating Roche lobe filling giant 
(Kcnyon et al. 1991). In T CrB, however, the observed 
vsuvi is by ~ 20 — 30% lower than expected for such 
a giant (Belczynski & Mikolajewska 1998| ). A very low 
wsini = 4.5 ± 2kms~ 1 has been also measured for the 



red giant in BF Cyg (Fekel et al. 2001). The problem is 
further complicated by the fact that the ellipsoidal vari- 
ations can be hardly visible only in quiescent visual light 
curves, whereas systematic searches for such changes in 
the red and near-IR range where the cool giant domin ates 
the continuum light are missing (Mikolajewska [2001 ). In 



Note that this estimate is based only on the shape of light 
curve and it is independent of the model adopted in their paper. 

3 The ratio of the observed red giant's radius to the mean 
tidal radius, R t is a better measure of the Roche lobe filling 
factor than the ratio of the giant's radius to the distance of the 
inner Lagrange point from the giant's center adopted by S01. 
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fact, a secondary minimum is also visible in the quies- 
cent optical light curves of AR Pav (Fig.l of Skopal et 
al. 2000a). Observations in the red and near-IR range are, 
however, necessary to confirm whether this minimum is 
due to ellipsoidal changes of tidally distorted red giant. 

Based on near IR TiO band depths, Miirset & Schmid 
( |1999| ) estimated a spectral type M5 III for the giant, 
whereas Allen ( |1980D deduced M6 III from the 2.3-^m CO 
band. If the giant is similar to normal giants, then its effec- 
tive temperature should be between T e j / (M5 III) = 3355 

199S| ) 



and 
and 



K and T e// (M6IH) = 3240 K (Richichi et al. 
luminosity, L g ~ 1800 - 4700 L Q , for R g = 130 R 
R g = R t = 210 Rq, respectively. Comparing the position 
of the red giant in the HR diagram with the evolutionary 
tracks for RGB/AGB stars, we find the red giant's mass 
range of 2M Q < M < 5 Mq (e.g. Vassiliadis & Wood 



199J; Schaller et al. |1992) Bessel et al. |1989|) , consistent 



with our dynamical mass estimate. 

The radial velocity curve of the F-type absorption lines 
suggest that the shell spectrum is formed in the outer- 
most accretion disk and/or the region where the stream 
encounters the inner accretion disk or an extended en- 
velope around the hot component. The broad emission 
wings could be then formed in the inner accretion disk 
or the extended envelope around the hot source. TH74 
came to similar conclusion based on the radial velocity 
behavior and the fact that the Balmer continuum served 
as a strong background for the cF absorption lines. We 
can place some limit on the size of the inner disk radius 
assuming Keplerian rotation. Then we have 

i? d /R w 19 x (M h /M )(w d /lOOkms _1 ) _2 sin 2 i. 

In particular, the Ha line wings follow the hot compan- 
ion's orbit down to at least ±270 km s -1 (b 6 A; Fig. |J), 
which corresponds to a Keplerian disk radius of 2.6 and 

o o 

2.8 Rq, for i = 90 and 70 , respectively. These values are 
much lower than the radius of the eclipsed object, i?h, C ci ~ 
0.1 A ( 46 R(p), de rived from optical eclipses (Skopal et al. 
2000a; Andrews 1974). Unfortunately, our method does 



not allow to measure accurately the wing positions for 
b ^ 6 A because of the presence of strong central absorp- 
tion. We also note, that the projected Keplerian velocity at 
an outer disk rim with R^ ~ 46 Rq is ~ 128kms _1 , very 
close to the average peak-to-peak separation in Balmer 
lines, ~ 148 for Ha and ~ 117 km s -1 for H/3, respectively. 
However neither the two peaks nor the central absorption 
trace the hot component orbital motion (see also Sec. 3.3). 
Thus we do not believe that the observed double-peaked 
structure originates in an outer accretion disk. 

The maximum velocity in the Ha line profile, 
800 km s -1 , is consistent with a Keplerian velocity at a 
distance of ^ 0.3 Rq, which sets an upper limit for the 
radius of the hot component. Although this value is an 
order of magnitude larger than the expected radius of a 
1 Mq white dwarf, it falls in the range of the observed 
hot component radii in classical symbiotic systems (e.g. 
Miirset et al. 1991). This result also rules out the pres- 



The permanent presence of the strong central absorp- 
tion in H I Balmer emission lines, stationary with respect 
to the orbital motion, indicates that the nebula is bounded 
on all sides by significant amount of neutral material (see 
Sec. 3.3). On the other hand, the strong double-peaked 
[O in] lines require an extended low-density ionized re- 
gion. A plausible solution to this apparent inconsistency is 
a bipolar ionization structure of the nebula, with a lower 
density towards the poles than in the orbital plane. We 
note here that similar strong [Om] lines together with 
strong central absorption in H I emission lines were ob- 
served during late outburst phases and decline in CI Cyg 
and AX Per. Moreover, in AR Pav, the eclipses in H I and 
He II lines are narrow with well-defined eclipse contacts 
(TH74; Quiroga et al., in preparation) likewise the eclipses 
in CI Cyg and AX Per dur ing the outbur st and its decline 
(Mikolajewska & Kenyon 1992a , 1992b ) which indicates 
that they may arise in a non-spherical, perhaps bipolar, 
nebula or flow. Such a geometry is also consistent with 
the re cent r adio observations of CI Cyg (Mikolajewska & 
Ivison 2001 ). These bipolar structures could be associated 
with the presence of an accretion disk in AR Pav and re- 
lated symbiotic systems at least during their active phases. 

Even if the giant does not fill its Roche lobe, its wind 
is likely focused towards the secondary and/or towards 



the orbital plane (Mastrodemos & Mor ris |1998{ , [1999 
Gawryszczak, Mikolajewska & Rozyczka |2002 ). In par- 



ticular, 3-D hydrodynamic model calculations show that, 
in a system with orbital parameters similar to those of 
AR Pav, gravitational interaction of the cool giant's with 
the secondary all alone can produce an equatorial to po- 
lar density contrast as large as 100-1000 giving rise to a 
bipolar geometry of the circumstellar nebula even in the 
absence of an accretion disk. 

A bipolar geometry of the nebula of AR Pav can also 
account for the orbital phase-dependent changes of the 
[Om] and [Nil] emission lines. In particular, such varia- 
tions can be expected if these lines are formed in the hot 
component wind confined by the asymmetric dense wind 
from the red giant o r in t he wind-wind interaction zone 
(e.g. Walder & Follini 200C). The radial velocity curves for 
the [O in] and [N n] emission lines, respectively, are not in 
phase, presumably because the two ions probe different 
regions, in particular with different density (note that the 
forbidden line emission is the most effective if the electron 
density is near to the respective critical density). 

4. Concluding remarks 

We have obtained the circular orbits for both components 
of AR Pav from direct radial velocity measurements of 
the M-giant and the hot component. In particular, we 
have found that the hot component orbital motion can 
be successfully traced by radial velocities derived from 
the broad emission wings of Hi Balmer lines and Hen 
lines. Similar result has been recently obtained by Ikeda 



ence of a main-sequence accretor in AR Pav. 



& Tamura ( |200C| ; also Fekel et al. |2001| ) for V1329 Cyg, 
in which the Hi and Hell broad wings appear to fol- 
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low the hot component orbit. With both radial velocity 
curves, we have obtained the orbital parameters and dy- 
namical masses for each component. The masses obtained 
are M g sin 3 i = 2.5 ± 0.6 M Q and M h sin 3 i = 1.0 ± 0.2 M Q 
for the red giant and the hot component respectively. 

Our M-giant orbit is in good agreement with that 
given by SOI. We believe that the somewhat higher semi- 
amplitude derived from combined our and SOI velocity 
data results from a better coverage of the final radial veloc- 
ity curve near its minimum and maximum. Our improved 
mass function, /(M g ) = 0.079 Mq, combined with the red 
giant's mass, M g = 2.0 Mq, derived by S01 from evolu- 
tionary tracks yields Mh = 0.87 Mq, in agreement with 
our dynamical mass estimate. As in S01, our final com- 
ponent masses are very different from earlier published 
values (TH74; Skopal et al. 2000a ). In particular, our up- 
per limit for the hot component mass, Mh 1.2 Mq, is 
much lower than the 2.5 Mq and 4.5 Mq values found by 
TH74 and Skopal et al., respectively, and it does not re- 
quire the presence of a main-sequence companion to the 
M giant. Moreover, our analysis (Sec. 3.4) practically rules 
out such possibility. Instead, we find that the hot compo- 
nent, like in most symbiotic stars, resembles central stars 
of planetary nebulae. However, we must stress that the 
hot component of AR Pav seems to be one of the most 
massive among thusfar studied symbiotic systems. 

AR Pav is also one of the most active symbiotic sys- 
tems. In particular, its optical light curve in addition 
to eclipses shows strong cycle-to-cycle variability (e.g. 
Andrews |1974| ; Bruch et al. [1994| ; Skopal et al. |2000aD 
which is related to the presence and strength of the 
cF absorption spectrum and a blue continuum (TH74; 
Hutchings ct al. 1983). Such a blue A/F-type shell spec- 
trum is typical for symbiotic systems during outbursts (ac- 
tive phases), in particular for those with multiple Z And- 
type eruptions, and it seems to be assoc iated with th e 
hot component (Mikolajewska & Kenyon 1992a . 1992h| ). 
Usually, the blue shell component disappears during a few 
years after the visual maximum, and a typical symbiotic 
system spends more time in qu iescence than in activity 
(Mikolajewska & Kenyon 1992a ). In the case of AR Pav 
however the F-type component remained visible, although 
weakened following the decline of the optical brightness, 
throughout all the period 1953-1973 covered by the ob- 
servations of TH74, as well as during our observations in 
1990-2000. Moreover, a warm blue continuum was also 



present in 1980-1982 (e.g. Hutchings ct al. 1983 ); unfortu- 
nately there are not published optical spectra which could 
confirm the presence of any blue absorption system. Thus, 
the hot component of AR Pav, unlike the hot components 
in typical symbiotic systems, seems to remain active most 
of time. There is also rather no doubt that this high ac- 
tivity is somehow related to variable accretion processes. 

Finally, AR Pav shares many properties with W Ser 
stars (or hyperactive Algols). In particular, both seem 
to have the accreting star embedded in an optically 
thick shell with an A/F-type spectrum, and strong emis- 
sion from circumstellar plasma. In the case of W Ser 



stars, there is strong evidence that the shell spectrum is 
formed mainly in a geometrically and optically thick ac- 
cretion disk seen nearly edge-on and in a gas stream (e.g. 
Andersen et al. |1988 and references therein). Kenyon & 
Webbink fll984| ) propose the same interpretation for the 
A/F- type shell spectrum observed in AR Pav, CI Cyg and 
a few other symbiotic systems during their active phases. 
Our analysis of radial velocity patterns in the F-type shell 
absorption as well as various emission lines in AR Pav 
supports such interpretation. The main difference between 
these symbiotic systems and W Ser stars is that the He n 
recombination lines, such as A 1640 A and A 4686 A are 
absent in W Ser stars whereas they are strong in AR Pav, 
CI Cyg and other related symbiotic stars (except during 
the maximum of the strongest optical outbursts). This 
seems to be related to different nature of the accretor hid- 
den inside the disk, which in W Ser stars seems to be a 
main-sequence B-type star whereas in AR Pav and other 
symbiotics the accretor is less massive, Mh 1 Mq , much 
hotter, and perhaps more compact. Both AR Pav and W 
Ser stars also have very complex and variable light curves, 
which can be accounted for by instabilities of a luminous 
disk component, and/or variations in mass transfer and 
circumstellar gas in the system. A complex circumstellar 
envelope is also suggested by variable polarization in both 
AR Pav and W Ser stars (Brandi et al. 2000] and Plavec 
1980, respectively). 

Although it is very premature to claim that all sym- 
biotic systems with multiple outburst activity do have 
tidally distorted giants, and - at least during active phase 
- accretion disks, whereas the non-eruptive systems do 
not, the former seems be the case for CI Cyg, YY Her, 
and a few other active systems. The latter is also true 
in the case of two non-eruptive systems, V443 Her and 
RW Hya, which do not show any evidence for ellipsoidal 
variability in their near-IR light curves, and their near-IR 
magnitudes show only small fluctu ations with the obser- 
vational errors (Mikolajewska et al. 2002 ). 

AR Pav certainly deserves further studies. For exam- 
ple, near infrared monitoring can provide important infor- 
mation about the process of mass transfer and accretion 
by confirming or excluding a tidally distorted giant. The 
structure of the nebula and the hot continuum source may 
be reconstructed and tested during the eclipses. AR Pav 
has not yet been detected in radio range, and the available 
upper limit Fe cm Si 14 mJy (Belczynski et al. 2000 ) is very 
far from the present detection limits. Multifrcquency ra- 
dio observations would give us information about the mass 
loss and the geometry of ionized region. We also desper- 
ately need good models of the structure and spectra of 
accretion disks and boundary layers for symbiotic stars. 
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